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We investigate Landau damping of Bogoliubov excitations in a dilute Bose gas moving in an 
optical lattice at finite temperatures. Using a ID tight-binding model, we explicitly obtain the 
Landau damping rate, the sign of which determines the stability of the condensate. We find that 
the sign changes at a certain condensate velocity, which is exactly the same as the critical velocity 
determined by the Landau criterion of superfluidity. This coincidence of the critical velocities 
reveals the microscopic mechanism of the Landau instability. This instability mechanism shows 
that a thermal cloud plays a crucial role in breakdown of superfluids, since the thermal cloud is a 
vital source of Landau damping. We also examine the possibility of simultaneous disappearance of 
all damping processes. 
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I. INTRODUCTION 

Superfluidity, one of the most startling phenomena 
in atomic and condensed matter physics, has fascinated 
both experimentalists and theorists since its first dis- 
covery in liquid Helium |l|, |^. Unfortunately, the 
inter-atomic interaction of liquid Helium was too strong 
to investigate the microscopic properties of superfluid- 
ity. However, the revolutionary realization of the Bose- 
Einstein condensation (BEC) in alkali gases [^|^|^ pro- 
vided us new fertile grounds for microscopic investiga- 
tions. The diluteness of the gases allows us to treat 
many problems, including the static and dynamical prop- 
erties of the BEC, in a microscopic way |^. In fact, 
a number of captivating properties of superfluidity have 
been revealed in dilute Bose gases such as quantized vor- 
tices H H and the Josephson effect Ha More- 
over, with the recent technological advancement in opti- 
cal laser physics, several remarkable experiments on the 
breakdown of superfluidity in optical lattices have been 
reported El Elii 15, 16J. 

Regarding the breakdown of superfluidity, it has been 
known that a superfluid current has its unique critical 
velocity, above which the current suffers friction, leading 
to the decay of superfluidity P, 0, ■ Landau pointed 
out that this decay is due to the appearance of negative 
excitation energy, inducing the energetic instability of the 
superfluid state [13 ■ Applying Landau's argument to the 
uniform system of a dilute Bose gas, the critical velocityis 
known to coincide with the Bogoliubov sound velocity [g ■ 
Although Landau originary assumed the homogeneity of 
the system, his idea of the energetic instability, which is 
known as Landau instahility, is also applicable to non- 
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uniform systems such as the BEC in optical lattices [23, 

El. 

In addition to the Landau instability, another kind 
of instabilit y i s passionately studied in recent years 
|l4i I15I I20I I21I |23 | . It is called dynamical instability, 
which is theoretically related to the appearance of the 
imaginary component of excitation energy (2Qi . i21.] , indi- 
cating the exponential increase of collective oscillations 
of the condensate. Thus the dynamical instability is ex- 
pected to be observed even at T = 0, while the Landau 
instability is not likely to be observed at T = 0, since 
it requires the thermal (non-condensate) component that 
receives the energy emitted by the condensate during the 
breakdown process. 

These two kinds of instabilities were experimentally 
studied for Bose gases in moving optical lattices in 
Refs. Ellli. Their results show that the energetic insta- 
bility takes place only at finite temperatures with a suffi- 
cient amount of thermally excited non-condensate atoms. 
This experimental result strongly suggests that thermal 
clouds are crucial to induce the energetic instability of 
the superfluid flow. 

In this paper we theoretically investigate the effect of 
the thermal clouds on the stability of the superfluid cur- 
rent. It should be noted here that this type of issue can- 
not be studied by the usual Gross-Pitaevskii equation 
containing no thermal component. 

We restrict ourselves to the collisionless regime. In a 
low density system at a low temperature, the mean free 
path of the excitations is known to become comparable 
with the size of the system itself and collisions between 
atoms play a minor role _23.] . 

Recently, Konabe and Nikuni j23| studied the instabil- 
ity of the condensate with a superfluid current in a ID 
optical lattice using a dissipative Gross-Pitaevskii equa- 
tion ^3 in the collisionless regime. They showed that 
the collisional damping of a collective mode due to the 



collisions between the condensate and the thermal com- 
ponent leads to the instability at a critical quasimomen- 
tum of the condensate, which exactly coincides with the 
Landau criterion. 

In contrast, in this paper we focus on the role of Lan- 
dau damping process in the instability of the conden- 
sate. Landau damping arises due to a coupling between 
a collective mode of the condensate and a thermal ex- 
citation. This damping process is known to be a dom- 
inant damping process in the coUisionless regime in the 
usual trapped Bose gas. In fact, Williams and Griffin 
[2fij | showed that the Landau damping rate is larger than 
the collisional damping rate, which is studied in Ref. ^^^I. 
Landau dam pin g has been studied by many theorists in 
uniform [23, EL trapped |22|, (p |^ ^1 l^d, and 
optical lattice |3J] systems. While all of these researches 
focused on current-free systems, in this present paper we 
study the Landau damping in a current-carrying dilute 
Bose gas and discuss the instability of the Bose-Einstein 
condensate induced by the Landau damping process. 

This paper is organized as follows. In Sec. El we de- 
velop the finite-temperature mean-field theory and derive 
the expression of the Landau damping rate for a current- 
carrying Bose gas. In Sec. IIII Al we develop the mean- 
field theory in a moving optical lattice within the tight- 
binding approximation. 



In Sec. IIII Bl we obtain an ex- 
pression for the damping rate for a Bose gas in a moving 
optical lattice, and we also examine the disappearance of 
all the Landau damping processes. In Sec. IIVI we ana- 
lyze the derived damping rate and discuss the stability 
of the condensate induced by the Landau damping pro- 
cess, where we find that the Landau damping process is 
vital for the Landau instability. In Sec.|Vl we discuss the 
similarity and difference between the Landau instability 
and the dynamical instability. 



II. FINITE-TEMPERATURE MEAN-FIELD 
THEORY AND LANDAU DAMPING 



In this section, we develop a finite temperature mean- 
field theory [s^ and introduce the Landau damping rate 
for a current-carrying Bose gas. The ground canonical 
Hamiltonian is given by 



K = 



H - fiN 

dr-0t(r,t)(^- 



2to 



where t4xt is an external potential, g 



a is the 



s-wave scattering length, and m is the atomic mass. The 
Heisenberg equation for the Bose field operator ?/'(r, i), 
which satisfies the Bose commutation relation, is then 



ih 



dt 



2m 



+ Vc^t - M V-Cr, t) 



+ 5^^(r,t)V;(r,i)V-(r,i). 



(2) 



Separating out the condensate part from the non- 
condensate part in the usual manner, the field operator 
is expressed as 



^(r,t) = $o(r)+V^(r,i), 



(3) 



where $o(r) — (V'(r,i)) is the time- independent, spa- 
tially inhomogeneous order parameter. Substituting the 
ansatz (O to Eq. (|2J), and using the mean field approxi- 
mation, we obtain 



2m 



M + Kxt + gino + 2n) 



$o(r) 



+ 5m$*(r) =0, 



(4) 



where we introduced the condensate, the non-condensate, 
and the off-diagonal non-condensate densities as 



"o(r) = |$o|', 
n(r) = (^t(r)^(r)), 

m(r) = {ip{r)^{r)), 



(5) 
(6) 
(7) 



respectively. Now we employ the Popov ansatz |3g: 

m(r) = 0, (8) 

which is known as a good approximation of a Bose gas at 
a finite temperature j23l |. In the Popov approximation, 
the equation for the condensate order parameter $o(r) 
(Eq. Q) is reduced to 



2m 



Ai + Kxt + g{no + 2fi) 



$o(r) = 0. (9) 



In a similar way, the equation for the non-condensate 
field operator ip{v,t) is found to be 



h Vc^t - M 

Zm 

+ 2gn{r)ip{r, t) 
+ 2g<^l{r)i,Hr,t). 



(10) 



The equation for '0t(r,t) is the Hermitian conjugate of 
Eq. (|f This pair of equations is solved by expressing 
ip{r, t) in the form 



-iEjt 



/'^-f^WaJe'^^H , (11) 



2 



where ctj and are the annihilation and the creation op- 
erators satisfying the Bose commutation relation. Substi- 
tuting Eq. l(TT|l and its Hermitian conjugate into Eq. (fllUl . 
we obtain the equation for the quasi-particle amplitudes 



and Vj: [Bogoliuhov-de Gennes equation) 



L 



-L 



where L is defined as 



L = — h Vcxt - + 2.gn(r). 

Zm 



(12) 



(13) 



The local density of the excitations h can be expressed 
in terms of the quasi-particle amplitudes as 



n{r) 



E 



3 V 



(14) 



where fj is the Bose distribution function of the jth ex- 
citation 



al3Ei 



1 



(15) 



Now we are ready to introduce the Landau damp- 
ing rate. Consider a low-lying collective mode (quasi- 
momentum p) of the condensate (quasi-momentum k) 
with energy E^^p^ surrounded by a static thermal cloud. 
The thermal cloud can absorb (or emit) quanta of collec- 
tive modes of the condensate, leading a damping of the 
collective oscillation. Here we focus on the interaction be- 
tween a collective mode of the condensate and a quasipar- 
ticle of the thermal cloud; the collective mode p is coupled 
with the thermal excitation j, and transformed into the 
thermal excitation j, or vice versa. This process is known 
as the Landau damping process and it has been studied 
for current-free cases 

(There is another damping process known as Beliaev 
damping process; however, because its contribution is 
much smaller than the Landau damping at finite tem- 
peratures [23, we do not consider the Beliaev damping 
here.) 

The Landau damping rate is calculated by perturba- 
tion theory [s^l and given by 



III. LANDAU DAMPING IN A DILUTE BOSE 
GAS MOVING IN AN OPTICAL LATTICE 

A. Finite-temperature tight-binding theory 

Li this subsection, we develop a finite-temperature Bo- 
goliubov theory for a Bose gas moving in a ID optical lat- 
tice within a tight-binding approximation. We consider a 
Bose gas in a cigar-shaped magnetic trap combined with 
a periodic optical lattice. The radial trapping frequency 
is assumed to be large enough to ignore the motion in 
this direction. In addition, we assume that the longitu- 
dinal trapping frequency is negligibly small so that we 
ignore the effect of the magnetic trap in the longitudinal 
direction. As performed in recent experiments [3 j 
we consider the optical lattice moving in the longitudinal 
direction at a constant velocity —v. In order to properly 
apply the argument developed in Sec. ^ we employ the 
reference system which moves along with the optical lat- 
tice at the constant velocity —v. In this reference system, 
the superfluid part of the Bose gas flows with the velocity 
-y, while the lattice potential is at rest and expressed as 



~d) 



where 



2m(P 



(18) 



(19) 



is the photon recoil energy, s is the dimensionless pa- 
rameter, and d is the lattice distance. We assume that 
iksT/Eji ^ s so that the energy gap between the first 
and second excitation band is much larger than the tem- 
perature, and thus the first band is thermally occupied. 
We also assume that the optical lattice well is deep 
enough so that the wave functions are well localized on 
each lattice cite, allowing for using the tight-binding ap- 
proximation [34 .37. .38J : 



V'(z,i)=^[^o£ + <^Ki)]/(2"^d), 



(20) 



Ik 



,p = 4V ^ |A,,,f (/, - Jj)5{Ek,p +E,- Ej). (16) 



^,3 



Since we consider the current-carrying case, Ai,j is 
a slight gen eralization of the expression presented in 
Refs. mil: 



A,. 



dr 



cf>o{Uj,(uiM* -I- ViVj) + VpUiU*} 



+ ^o{UpiUiU* + ViV*) + UpViU*} 



(17) 



Using these general expressions, we derive the damp- 
ing rate for a Bose gas moving in an optical lattice in 
Sec. lnTBl 



where f{z — Id) is a localized function at the ^th site, 
which satisfies the orthonornial relation: 



dzf{z - id)f{z - md) = (5^™, 



(21) 



where Si^m is the Kronecker's delta. With f{z — it), (j)og 
and (pi (t) form the condensate order parameter and non- 
condensate field operator at ^th site respectively. 

The following formalism is an analogue of the argu- 
ment developed in Sec. |n] Substituting Eq. into 
Eq. ((2Jl, and using the Popov ansatz (jHJ, we obtain the 
discrete Gross-Pitaevskii equation including the mean- 
field of non-condensate parts 
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FIG. 1: Condensate fraction no/n in an optical lattice as a 
function of the condensate quasimomentum fc/gs at different 
temperatures r\ — ksT/J, where gs = vr/d is the Bragg vector 
and U 1 3 is fixed to 1. 



-M0o^ = 0, (22) 

where condensate densities noe, non-condensate densities 
rif, the hoping matrix elements J, and the on-site inter- 
action energy U are defined as 



At/ 



hi{z) = {(j)\{z)<j)i{z)) , 



J = - dzf, 



t+ii 



U^g j dz\f, 



(23) 
(24) 

(25) 
(26) 



respectively. Here we introduced g — 2h^a/ma\ as the 
coupling constant modified for our eff'ective ID system 
using the radial length aj_ 

Eq. (|22|l admits the solution of the Bloch form (We set 
h — \ from now on.) 



iked 



(27) 



where k is the quasimomentum of the condensate. Sub- 
stituting Eq. (|27|l into Eq. (|22|l . one can write the chem- 
ical potential /i and the order parameter $o 



fi = -2Jcos(fcd) + U [no(T) + 2no{T)] , (28) 



and 



e 



where Uq is the mean number of the condensate particle 
at a lattice cite. Normal mode solutions are found by 
expanding (j)e{t) as 



Mt) = E 



ikld 



(30) 



Substituting Eq. H3U|I and its Hermitian conjugate into 
Eq. H22|) , we obtain a pair of equations 



2Jsin(/cd) sin(pd) -I- 4Jcos(fcd) sin^(pd/2) 



k,p-U(l}oeVk,p= Ek,pUk^p, (31) 



2Jsm{kd) sm{pd) — cos{kd) sm^{pd/2) 



Uno 



vip + U^*oiulp = Ek,pilp, (32) 



where the excitation energy E^^p is 



Ek^p = 2Jsm{kd) siTi{pd) + Ek,p, 



(33) 



with 



El; 



k,p 



4j2cos2(fcd)sin4 ' 



+ 2JU TiQ cos(fcc?) sin 



2 I Pd 



(34) 



The Bogoliubov amplitudes {u, v), which satisfy the nor- 
malization condition 



dz[ \ Up - \ Vg p 



^p,qi 



(35) 



are expressed in Bloch forms as 



Uk.'j 



.i(P+k)edr( _ 



fiz - id) (36) 



and 



Vk,p{^) = ^ E ^Ip-^^^ " ^'^^ 



= ^T4^,^e'(f-'=)"/(z-^d). (37) 

In these expressions, / is the number of lattice cites, while 
Uk,p,Vk,p are given by 



(29) Uk^p 



\ 



4Jcos(fcd) sm^{^) + Una 



E, 



H (38) 



4 



\ 



4Jcos(fcd) sin^(^ 



[/no 



Ek.r. 



These satisfy the normahzation condition 



Uk,p 



1, 



(39) 



(40) 



j/ Ek^p is a real function. We note that | Uk,p P ^ | 
Vk.p P= if Ek,p contains an imaginary component. This 
problem was first pointed out by Niu and Wu [23|- We 
confirm that it is stiU vahd in our tight-binding model. 
We will discuss this important issue again in Sec. El re- 
garding the instability of condensates. 

The number of the condensate at a lattice site, no, is 
determined by the relation n = no -I- n, where n is given 
in Eq. H14|) . Note that no implicitly depends both on 
the temperature T and the condensate quasimonientum 
k. Fixing parameters U/ J — 1, n — 2 and / = 250, 
we solved the equation n = uq + n self-consistently. In 
Fig. ^ the condensate fraction no /n at different temper- 
atures rj = k^T/J is shown as a function of k. Due to 
the fatal limitation of the Bogoliubov theory, one cannot 
determine no at fc > A:c, where kc is the critical quasimo- 
mentum determined by Landau criterion of superfluidity 
[2(il |. since the low-energy excitation becomes negative. 
We will give the explicit expression of kc in Sec. lIVBl af- 
ter showing in Sec. lIV Al that kc turns out to be the same 
as the critical quasimomentum determined from the Lan- 
dau damping rate. 



B. Landau damping in a moving ID optical lattice 

Now we calculate the Landau damping rate for a Bose 
gas in a moving ID optical lattice using the tight-binding 
approximation introduced in Sec. IIII Al First, we cal- 
culate the matrix element Ai_j in Eq. (|17|l . By using 
Eqs. (113), ^ and we obtain, 



A, 




'p+q,.q,+G„ 



-t- VpU.Uj + Vp{U,U, + V,Vj) + UpV.Uj 



(41) 



where G„ = 2mr/d (n is an integer) is the reciprocal 
lattice vector corresponding to the Umklapp scattering 
process when n is not zero. We see that Eqs. and 
(|41|l involves the momentum and energy conserving delta 
functions. This means that, in order for the Landau 
damping to occur, there must exist qi and qj that satisfy 
the energy and momentum conservation 34] . Combining 
the two conservations, we require that there exists q that 
satisfies 



Ek 



k.p 



E, 



k,q 



E, 



k,p+q+G„ 



(42) 



CO 




k/q. 



FIG. 2: Quasimomentum p satisfying Ep + Ep+p — Ep with 
p <^ p. Both p and k are normalized by the Bragg vector 
QB = 7r/d, where d is a lattice spacing, a = Uno/J is set to 
1. 



For small quasimomentum <C 1, -Bfe,p in Eq. H33() has 
a linear dispersion in terms of p; namely 



Ek,p = [-2Jsin{kd) + y/2JUnoCos{kd)] \p\, (43) 

where we set the fc's direction as positive and considered 
the energy dispersion with p < 0, since only this type of 
excitation is responsible for the energetic instability, as 
we show in Sec. IIVI 

Assuming that p is small enough that Ek,p can be ex- 
pressed as Eq. (|43|) , one can reduce required the requiring 
condition Eq. 142|l reduces to 



[2Jsin{kd) - y/2JUno cos(fcd)] 



dEp 
d{pd) ' 



(44) 



where we denoted p = q ^ Gn- Note that the Umk- 
lapp scattering is a mere superficial problem for us here, 
since p only appears in the periodic energy expressions 
in Eqs. ^ and 

Fig. m shows p in Eq. 1)44(1 as a function of k with 
a = C/no/J = 1. It is seen that the Umklapp scattering 
appears at a finit e- fc, while it cannot be seen when k = Q 
as shown in Ref. |3J| . 

The condition for p to exist, determined by Eq. H44|l . 
is shown in Fig. |2| At fc = 0, all the Landau damping 
processes disappears when a > 6, which agrees with the 
result derived in Ref. With increasing fc, we find 

that the region where the Landau damping process dis- 
appears prominently decreases. This disappearance can 
be confirmed in future experiments. 

By integrating Eq. (|16|l using the delta-function for the 
energy conservation, we explicitly obtain the expression 
of the Landau damping rate for a Bose gas moving in a 
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0.5 




FIG. 3: Damping possibility diagram. Damping processes 
are possible to exist only in the shaded area. Dashed line 
indicates the critical quasimomentum of the condensate given 
by Eq. (HSl . 

ID optical lattice: 

Ek p UPd? JUno cos(fcd) 
^ \d'E,,,/d^p\ 4sinh^(^) 

"8Jcos(fcd) sin^r^) + Uno 
X = 

-1 2 

J(C/no) cos(fco?)sin(p(i) ^^^^ 
El py/ JUno cos(fcd) 

IV. LANDAU DAMPING AND LANDAU 
INSTABILITY 

A. Landau damping and Landau instability 

In this section, we show that the Landau damping pro- 
cess is closely related to the breakdown of superfluidity 
due to the energetic instability (or Landau instability). 
As explicitly shown in Ref. |23, the amplitude of the 
collective mode S^k,p decays (or grows) as 

S^k,p (X e-T^-*, (46) 

where jk.p is the Landau damping rate discussed in this 
paper. This relation indicates that the condensate with 
a superfluid current k is stable as long as ^k,p is positive 
for any collective mode momentum p. In fact, in a uni- 
formly resting (fc = 0) condensate, the Landau damping 
rate is always positive. Thus an induced collective mode 
decays exponentially in time [23l l30l | . This means that 
the thermally populated cloud stabilizes the condensate 
part. However, as shown below, it is possible that jk.p 
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FIG. 4: Landau damping rate yk,-o.oiqB ™ the ID optical 
lattice with U/J = 1 as a function of a condensate quasimo- 
mentum fc < fee at different temperatures — kBT/J. j and 
k are normalized to the hopping matrix element J and the 
Bragg vector qb respectively. 

becomes negative for a moving condensate with finite- 
k. Negative "fk,p indicates the appearance of the inverse 
process of the usual Landau damping process, whereby 
the thermal clouds emit excitation quanta. The increase 
of the amplitude of collective modes severely destabilizes 
the condensate state, leading to the breakdown of the 
superfluidity. 

We now show that the Landau damping rate in 
Eq. (|45|l changes its sign at a critical quasimomentum. 
We find that the sign of ^k,p fully depends on the sign of 
Ek,p, i-e. 

sgn(7fc,p) = sgn{Ek,p). (47) 

This means that the critical quasimomentum of the Lan- 
dau damping rate coincides with the critical quasimo- 
mentum of the corresponding collective mode determined 
by the Landau criterion, which is the Bogoliubov sound 
velocity when we consider a dilute Bose gas in a uniform 
system. Landau originally argued that the superfluid 
decays by emitting phonon, whose energy level is lower 
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than the superfluid state, which is indeed possible when 
the superfluid possess a velocit y l arger than the criti- 
cal velocity (Landau criterion) Here we find that 
the critical velocity obtained from the Landau damping 
rate coincides with the velocity derived from the Landau 
criterion, meaning that the Landau damping process is 
the microscopic mechanism of the macroscopic Landau 
instability. 

This close relation between the Landau damping and 
the Landau instability (Eq. (|47|l ) can also be derived from 
the general expression of the Landau damping rate in 
Eq. 11611 . The sign of jk,p in Eq. is determined by 
the difference fi — fj between the Bose distribution func- 
tions of static thermal clouds. Combining this with the 
energy conserving delta-function, we find that "/k,p be- 
comes negative when Ek^p becomes negative (Eq. (071). 
This means that the Landau damping process is generally 
the source of the Landau instability within the approxi- 
mation of static thermal clouds. 

Fig. 01 shows the Landau damping rate at k < kc at 
different temperatures 77 — ksT / J. The damping rate 
prominently decreases with increasing fc, meaning that 
the stabilization mechanism of the condensate by the sur- 
rounding thermal cloud (usual Landau damping) works 
less effectively at a larger k < kc- 
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FIG. 5: Stability diagram of the BEG in the optical lattice 
with a = 0.5. The shaded area is where the excitation en- 
ergy i5fc,p is negative thus the system is energetically unstable. 
In this area, the Landau damping rate 7fc,p is also negative 
(Eq. MTU ) and the inverse process of the usual Landau damp- 
ing is dominant. Both k and p are normalized by qs = 7r/d, 
where d is a lattice spacing. 



Although we do not show the growth rate 7fc^p of an 
unstable mode at fc > /cc due to the fact that the con- 
densate density at the region cannot be calculated in our 
Bogoliubov theory, we estimate its behaviour as follows. 
As suggested by the experiment in Ref. 0|, assuming 
that the condensate density remains constant for k > kc 
to its minimum value shown in Fig. ^ we confirmed that 
"fk,p prominently decreases with increasing k > kc and 
diverges to —00 at the edge k = Qb/^. Moreover, the 
decreasing rate of ^k,p is confirmed to be the increasing 
function of the temperature within the assumption. This 
indeed means that the amount of thermal clouds is cru- 
cially important in the breakdown process of superfiuids. 

To summarize, we propose the following general sce- 
nario for the breakdown of superfluidity due to the Lan- 
dau damping process with thermally populated excita- 
tions: As long as fc < fee, the Landau damping rate jk.p 
is positive, and thus the superfluid state is stable. How- 
ever, once k exceeds kc, the amplitude of collective modes 
increase exponentially in time by the inverse of process of 
the usual Landau damping process. This critical process 
should lead the breakdown of the superfluidity. On the 
other hand, if there is no thermal cloud, neither the Lan- 
dau damping nor the coUisional damping 24] can emerge. 
Thus the superfluid does not breakdown. 



B. Critical current 

We now evaluate the specific critical quasimomentum 
kc of the condensate, where the superfiuid starts to decay 
due to the inverse process of the usual Landau damping 
for a dilute Bose gas moving in an optical lattice. The de- 
sired kc is obtained by requiring dpE^^ ,p ii^ Eq. . 
Then, kc is expressed as a function of the quasimomen- 
tum p of an excitation; however, the critical current of 
a superfiuid should be the minimum current with which 
the superfluid starts to decay. As seen in the Fig. |S1 the 
smaller the quasimomentum of an excitation, the smaller 
the quasimomentum of the condensate required for nega- 
tive excitation energy to appear. Thus the critical quasi- 
momentum of the condensate should be the quasimomen- 
tum with which the energy of the phonon corresponding 
to the longest wave length becomes zero (see Fig.O. At 
this critical current, which is nothing but the Landau 
critical current, the superfluid state becomes unstable if 
the system contains enough amounts of thermal excita- 
tions. 

Therefore, by equating Eq. (|44|l to zero, we obtain 



cos{kcd) 



: + + 16 



(48) 



In short, we theoretically proved that thermal exci- 
tations play a significant role in the breakdown of the 
superfiuid due to the energetic instability. 



where a — Uhq/J. We note that the Landau instability 
is found to take place only where the Landau damping 
process exists as seen in Fig.O This shows that the very 
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close relation between Landau instability of condensates 
and the Landau damping process. 

We note that our ID results can be tested in the defini- 
tive ID optical lattice system using the recent technique 
[T^ I39I l40l| . which realizes a much tighter (in the ra- 
dial direction) magnetic trap than used in Ref. Ex- 
periments should also be conducted in the coUisionless 
regime. 

V. LANDAU INSTABILITY AND DYNAMICAL 
INSTABILITY IN AN OPTICAL LATTICE 

In this paper, we showed that the Landau damping 
works as a source of Landau instability at a finite temper- 
ature. The thermal cloud around the superfluid tempts 
the condensate into exciting by accelerating an energet- 
ically irresistible process, i. e. the inverse process of the 
usual Landau damping. This explains the microscopic 
mechanism of the Landau instability at a finite temper- 
ature. 

There is another well known instability of the con- 
densate in an optical lattice, the dynamical instability 
[20I I21I . In this section, we briefly discuss how the 
two kinds of instabilities are understood in connection 
with our results. 



A. Similarities 

As described in Sec. IIVI a collective mode S^k,p of the 
condensate grows or decays in time as 

S^k,p oc e~*-^'=-f*-'''=.p*. (49) 

If 7fc,p becomes negative, the corresponding collective 
mode oscillation exponentially grows, leading to the 
breakdown of the superfluid. This is Landau instabil- 
ity discussed in this paper. On the other hand, it is 
also possible in some systems that Ek^p itself contains an 
imaginary part. Then likewise, the excitation exponen- 
tially grows and the superfluid severely decays. This is 
the dynamical instability, which was characterized by the 
exponential growth of collective modes in time in former 
literatures |20l l2lj| . However, as seen in our argument. 
Landau instability can also be related to the exponential 
growth of collective modes. 

B. Differences: collective modes and elementary 
excitations 

One obvious difference between Landau instability and 
the dynamical instability is the scales of the critical veloc- 
ity, at least for a Bose gas in an optical lattice discussed in 
this paper. As seen in Eq. (|33|l . Landau instability takes 
place at < A: < (7b/2 while the dynamical instability 
occurs at (7b/2 < fc, where Ek,p contains an imaginary 
component. 



However, we note here that the most significant dif- 
ference between Landau instability and the dynamical 
instability is weather it is caused by modulational collec- 
tive modes or quantum elementary excitations. We use 
the collective mode for the "classical" modulation of the 
condensate itself, while the elementary excitation is used 
to mean the "quantum" and thermal depletion of the 
non-condensate part 6] . 

As seen in Eq. H40() . the Bogoliubov amplitudes {U, V) 
satisfy the usual normalization condition as long as the 
corresponding energy Ek.p in Eq. is a real function. 
Then (U, V) can be sufficiently quantized and Ek^p can 
be regarded as the energy of an elementary excitation 
However, once E^^p contains an imaginary part, or 
equivalently E^^p becomes purely imaginary, {U, V) can- 
not satisfy the normalization condition, but instead (with 
some straightforward algebra) 

I Uk,p P - I Vk,p P= 0, (50) 

which is also obtained using an abstract model in 
Ref. poj . This means that {U,V) cannot be quantized, 
thus they should be regarded as the amplitudes of the 
collective modulation of the condensate. 

In summary. Landau instability is caused by the ther- 
mal elementary excitation while the dynamical instability 
is caused by the collective modulation of the condensate 
itself. Consequently, the Landau instability only occurs 
at finite temperatures while the dynamical instability can 
be induced even at T = 0. This agrees with our analy- 
sis on the Landau damping that there should be a ther- 
mal cloud in order to induce the Landau instability since 
the thermal cloud is also a vital source for the Landau 
damping process. As noted in Sec. this temperature 
dependence of the decay of superfluids is experimentally 
observed in Ref. [Tsf . which is consistent with our anal- 
ysis on the elementary excitations and collective modes. 
We hope that our work will stimulate further experiments 
on the breakdown of superfluids. 



VI. CONCLUSION 

In summary, we investigated the instability of a mov- 
ing condensate in an optical lattice surrounded by ther- 
mally excited quasiparticles at finite temperatures. We 
revealed that the thermal clouds destabilize the conden- 
sate by the inverse process of the usual Landau damping. 
Also, we explicitly calculated the Landau damping rate 
for a Bose gas at a finite temperature and showed the re- 
lation between Landau damping and Landau instability. 
We found a certain condition for energy scales that Lan- 
dau damping or Landau instability occurs. We finally 
noted the difference and similarity between the Landau 
instability and the dynamical instability in an optical lat- 
tice. 
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